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The scheme recently proposed in [M. Scala et al, Phys Rev Lett 111, 180403 (2013)], where a 
gravity-dependent phase shift is induced on the spin of a nitrogen-vacancy (NV) center in a trapped 
nanodiamond by the interaction between its magnetic moment and the quantized motion of the 
particle, provides a way to detect spatial quantum superpositions by means of spin measurements 
only. Here, the effect of unwanted coupling with other motional degrees of freedom is considered and 
we show that it does not affect the validity of the scheme. Both this coupling and the additional 
error source due to misalignment between the quantization axis of the NV center spin and the 
trapping axis are shown not to change the qualitative behavior of the system, so that a proof-of- 
principle experiment can be neatly performed. Our analysis, which shows that the scheme retains 
the important features of not requiring ground state cooling and of being resistant to thermal 
fluctuations, can be useful for the several schemes which have been proposed recently for testing 
macroscopic superpositions in trapped microsystems. 


I. INTRODUCTION 

Since its discovery more than one century ago, quan¬ 
tum mechanics has puzzled the scientific community with 
questions foundational. Indeed, while everybody agrees 
on the power of quantum mechanics for the description of 
microscopic systems, such as atoms and molecules, and 
on its applicative power, from electronics to the most 
recent developments in quantum information processing 
and computing [1], there is still a lot of debate about 
the transition from the microscopic to the macroscopic 
world, where experience shows that quantum mechanics 
does not seem to be valid and is substituted by classical 
physics [2]. As an example, the naive use of quantum 
mechanics for the description of the macroscopic world 
would lead to predictions contradicting our experience, 
such as in the well known case of Schrodinger’s cat [3]. 
While it is a relatively simple task to define what the mi¬ 
croscopic and the macroscopic worlds are, it is difficult 
to draw a line separating these worlds. Finding the bor¬ 
der between the two worlds is an important foundational 
issue and, so far, many solutions have been proposed to 
describe the transition from quantum to classical physics, 
from more widely accepted theories based on decoher¬ 
ence processes due to the interaction with an external 
environment to more debated theories, which pro¬ 

pose modifications of quantum mechanics but have not 
been tested so far, such as the spontaneous localization 
theories [Ms] 

In this framework, mesoscopic physics, i.e. the physics 
of systems which lie somewhere in between the micro¬ 
scopic and macroscopic worlds, plays an important role, 
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since one can play with parameters which are intuitively 
related to the transition from quantum to classical, such 
as the total mass or the total number of atoms involved 
in the dynamics of the systems. Examples of experiments 
with mesoscopic systems are given by double-slit inter¬ 
ference with very large molecules [9], the production of 
nonclassical states of light by means of optomechanical 
systems [TOHTS] and the study of coherence in trapped 
nanoparticles mHiH]. In all these experiments, the scien¬ 
tific community is producing quantum superpositions of 
states of larger and larger systems, and there has been a 
lot of work about quantifying the macroscopicity of such 
superpositions [IMB- The generally accepted idea is 
that, by taking more macroscopic regimes and by making 
the interaction with the environment weaker and weaker, 
we can finally reach a regime wherein alternative theo¬ 
ries, as opposed to orthodox quantum mechanics, will 
become experimentally testable [8]. 

Usually, the experiments to test such macroscopic su¬ 
perpositions can be quite involved, since in general they 
require the possibility of cooling the system, coupling it 
with cavity resonators, and also resolving the spatial ex¬ 
tension of the particles under study [22]. Moreover in 
general one will need an ensemble of identically prepared 
systems. In our previous work we proposed a quantum 
interference scheme for trapped nanoparticles which over¬ 
comes these requirements and would allow for the detec¬ 
tion of the quantum features of the motion of the particle 
by spin measurements only [23|. In the proposal, we con¬ 
sidered a conditional displacement induced by a magnetic 
field gradient on a trapped nanodiamond containing one 
nitrogen-vacancy (NV) center. Starting from a superpo¬ 
sition of two distinct states of the spin of the NV center 
and from a generic coherent state for the motion of center 
of the mass (CM) of the trapped bead, we showed that 
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the phase difference acquired by the different trajecto¬ 
ries of the harmonic motion in the presence of gravity 
can be completely transfered to the spin states, so that 
standard Ramsey interferometry is all we need to detect 
the gravity-induced phase difference. 

The analysis performed neglected some terms in the 
Hamiltonian making it effectively one-dimensional, in¬ 
stead of genuinely three-dimensional as it would be in real 
experiments with optical tweezers. The simplified model 
allowed us to catch the features of the scheme proposed 
in an easily understandable way. In this paper we will 
show that a perturbative approach can be used to prove 
that those additional terms can be neglected up to very 
high values of the relevant coupling constants, which im¬ 
plies that the confinement along x and y need not be as 
tight as one would naively think: this is due to the fact 
that the first order correction to the energy due to the 
additional terms are zero and so one has to go to second 
order to get corrections to the dynamics. We will indeed 
show that, for parameters describing correctly the dy¬ 
namics of a nanodiamond trapped in an optical tweezer, 
the fidelity between the perturbatively corrected state 
and the state predicted by the one-dimensional model is 
always larger than 99%. The same perturbative tech¬ 
niques can be used to treat the effect of misalignment 
between the quantization axis of the spin of the NV cen¬ 
ter, which essentially depends on the relative direction 
between the color center and the diamond lattice, and 
the trapping axis z. We will show that misalignment 
does not change the qualitative behaviour of the system 
dynamics and that we are still able to detect interference 
fringes showing the gravitationally-induced phase differ¬ 
ence between the spin states. Moreover, the most im¬ 
portant feature of the one-dimensional scheme, i.e., the 
robustness of the phase difference against thermal fluc¬ 
tuations, is not substantially affected by the additional 
terms in the Hamiltonian: therefore the more realistic 
scheme shows that no (or a very small amount of) cool¬ 
ing is required in order to get the interference fringes, so 
that the simplicity of the proposal remains. The paper 
is structured as follows. In Section H we review the one¬ 
dimensional scheme, showing how the quantized motion 
of the particle leaves phase-signatures on the states of 
the spin of the NV center, while the effect of the coupling 
with the motion along the x and y directions and the ef¬ 
fect of misalignment are presented in Sections HI and IV 
respectively. In Section V the results are discussed and 
some concluding remarks are given. 


II. GRAVITATIONAL-INDUCED PHASES ON 
SPIN STATES 

A. The system and its dynamics 

As shown in Figj^ the setup consists of a nano-scale 
diamond bead containing a single spin-1 NV center lev¬ 
itated by an optical tweezer in ultra-high vacuum. The 



FIG. 1. An optical trap holds a diamond bead with an NV 
center with both weakest confinement and spin quantization 
along the z axis. A magnetized sphere at zo produces spin- 
dependent shifts to the center of the harmonic well. An angle 
0 between the vertical and the z axis places the centers of the 
wells corresponding to the | + 1) and | — 1) states in different 
gravitational potentials. Starting with an arbitrary coher¬ 
ent state, the CM of the bead oscillates as different coherent 
states in the centre-shifted, spin dependent well (red solid and 
dashed line), accumulating a relative gravitational phase dif¬ 
ference due to the superpositions. At to = 21^juJz this phase 
can be read from Ramsey measurements on spin. The blue 
shaped zone shows a generic orientation of the NV center’ s 
axis z with respect to magnetic direction z 


motion of the bead is coupled to the S = 1 spin of the NV 
center by means of a static magnetic field gradient which 
can be generated by a magnetized sphere with a perma¬ 
nent dipole moment m = (0,0, m;^) oriented along the ^ 
direction. Defining a reference frame such that the cen¬ 
ters of the harmonic potential and the magnetized sphere 
be at (0,0,0) and (0,0, zq) respectively, we can expand 
the magnetic field of the sphere around (0, 0, 0), and get 

B^ = -Box, By =-Boy, B,, = +2Bo z, (1) 

where Bq = 3jiiomz/{4:7rzQ). Therefore the interaction 
between the spin of the NV center and the vibrational 
motion can be described by the Hamiltonian 


ifint = -A [2 5, (c+ct)-d— 5, ia+a^)-J‘^Sy (6+6^)], 

V Y 

( 2 ) 

where 


A = 


SyorrizZo 

4:7t\zo\^ 


QNV 


h 


2muJz^ 


(3) 


m being the mass of the bead, qnv the Lande factor 
of the NV center, (Ib the Bohr magneton and a, 6, c 
the annihilation operator of the oscillation in x, z re¬ 
spectively. In this Section we will neglect the interac¬ 
tion between the spin and the x and y directions in Eq. 

whose effect will be analyzed in the next Section, on 
the basis that ^ In support of this approx¬ 

imation, in Fig. we present experimental data mea¬ 
sured in our laboratory from a nanodiamond levitated in 
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moderate levels of vacuum of approximately 10 mB us¬ 
ing 200 mW of 1064 nm trapping power. We measure 
^ ~ ^ 0.18. Due to asymmetry in the laser focus 

[25] . oscillation frequencies in the radial directions are 
separated by approximately ci;/27r = 5 kHz. The lower 
axial z frequency arises from the smaller electric field gra¬ 
dient along the beam axis in comparison to x and y. The 


U)b4 nm Vanmim chamber Balanced 

laser nhotodiode 



a) 



a) /2ji (kHz) 

b) 


FIG. 2. (a) Experimental schematic of the optical dipole trap. 
A 1064 nm laser beam is tightly focused by a high numerical 
aperture (0.95) objective. The polarization of the trapping 
light can be rotated by a half-wave plate. Scattered light from 
levitated nanodiamonds is collected by a lens and sent to a 
balanced photodiode in an interferometric scheme described 
in reference m, providing a position dependent signal from 
the levitated nanodiamond, (b) Power spectral density (PSD) 
as a function of u at approximately 10 mB using 200 mW of 
trapping power. Fourier transforming the position dependent 
signal yields the PSD of the trapped nanodiamond. The axial 
z frequency has been scaled by a factor of 20 for clarity. 

X OT y frequency is revealed to the balanced photodiode 
by rotating the polarization of the trapping light, while 
2 : is measured as the nanodiamond traverses an inten¬ 
sity gradient along the beam axis, giving rise to a time 
varying scattering intensity. 

Finally, we add the free Hamiltonian of the bead and 
of the NV center and the Hamiltonian describing the in¬ 
teraction between the bead and the Earth’s gravitational 
field, i.e., mgz cos 0, where 0 is the angle between the 2 ; 
direction (see FigQ, so that the total Hamiltonian of the 
system is given by 

H = DSl + huj;,c'<c-2{XS^-AX){c + €"<), (4) 

with 

AX=-mqcos0\ —-—. (5) 


The Hamiltonian above represents a harmonic oscillator 
whose center depends on the eigenvalue of Sz- In its 
derivation we have also assumed that the Zeeman split¬ 
ting of I-1-1) and | —1) due to the zeroth order expansion 
of Bz is cancelled by addition of a uniform magnetic field 
along z. 

It can be shown that, starting at t = 0 from the state 
\f3) \sz), where \^) is a coherent state of the center of mass 
(CM) quantized motion of the bead in the harmonic well 
centered at 2 : = 0 and \sz) is an eigenstate of the operator 
Sz with eigenvalue Sz = +1, 0, —1, the system evolves at 
time t to the state \/3{t,Sz)) Is;^) where 

\P{t,Sz)) 

X -h u) , (6) 

and u = 2{sz A — AX)/huJz- It is worth noting that, at 
time to = 211 jujz^ the oscillator state returns to its origi¬ 
nal coherent state for any /3 and Sz- One can take ad¬ 
vantage of this feature to show that spin measurements 
at to will be unaffected by any thermal randomness in 
the initial motional state of the oscillator. 


B. Detecting the gravitational field by Ramsey 
interferometry 

What we want to do now is to prepare the spin in a 
superposition of states, so that when the center of mass 
of the bead undergoes a conditional displacement due to 
dynamics, the hybrid system evolves to a state which in¬ 
volves superpositions of |/3(t, Sz)) |+1) and |;d(t, Sz)) | —!)• 
We then measure the spin at a special time to at which 
the conditional displacements are again undone due to 
the natural dynamics so that the spin state becomes un¬ 
entangled from the vibrational one. In this way we can 
perform measurements on the spin only, to reveal the 
different phases acquired by the spin components, due to 
the evolution for the vibrational state. 

Preparing the system in the separable state \/3) \sz = 0) 
and applying a microwave pulse corresponding to the 
Hamiltonian Hmw = (1+1) (0| + l~l) (0| + h.c.), with 

O much larger than any other coupling constant and for 
a pulse duration tp = 7r/(2\/2r^), the spin state becomes 

|T(0)) = 1/3) which we will take as the ini¬ 

tial state for the interaction under the Hamiltonian 
After the interaction time t, the state is 

|^ft,+l))|+l) + |/3ft,-l))|-l) 

^/2 

which is the superposition we intend to evidence. From 
the expressions of |/3(t,±l)) in Eq.([^ we can see that 
separated coherent states are involved in the above su¬ 
perposition along with phases due to gravitational poten¬ 
tial, which will evidence the above superposition. Up to 
a global phase factor, the state after an oscillation period 
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to = 211 /ujz is 


l^(io)) = 1/3) 




I_ 


1 ) 


x/2 


with 


A^Grav 


16AAA 

h?ujz 


to. 


( 8 ) 

(9) 


The phase difference A^crav can be revealed by applying 
Hmw again. Indeed, after a time tp, the population of the 
spin state with 5';^ = 0 is: 


Po{t 


to Ptp) = cos^ 


f A0Gr 


( 10 ) 


which gives a direct connection between the value of the 
phase shift and spin population. As A^Grav oc g can 
never appear as a relative phase between spin states un¬ 
less spatially separated states of CM were involved in 
the superposition (|;d(t,+1)) |+1) + |/3(t, —1)) | —1)) /V2 
for 0 < t < to, the detection of A^Grav evidences such a 
superposition. 

Finally, let us consider what happens if the initial 
state is the product of a thermal motional state pth 
and an eigenstate of the spin operator Sz- In fact, 
one can exploit the fact that the results given above 
are independent of the amplitude (3 and that any ther¬ 
mal state pth of the motion can be written as pth = 
f d^/3Pth(^) I/d) (/d|, where Pth is the Glauber P repre¬ 
sentation for the thermal state, to show that, after the 
evolution over one oscillation period, the state of the sys¬ 
tem is again factorizable and that the phase difference 
accumulated by the spin states is not affected by the 
thermal motion. Basically, though a mixture of many 
Schrodinger cats |/3(t, -hi)) |+1) + |/d(t, —1)) | —1) is gen¬ 
erated for 0 < t < to, the interference between the com¬ 
ponents |;d(t,+1)) 1+1) and |/3(t, —1)) | —1) of the cat is 
independent of /3. This immunity of the interference to 
thermal states hinges on the mass being trapped in a har¬ 
monic potential. We assume that anharmonic effects of 
the trapping potential will be avoided by feedback cool¬ 
ing of our oscillator to mK temperatures [24l [25] . 


III. THE EFFECT OF THE TERMS 
NEGLECTED: PERTURBATIVE ANALYSIS 

In the previous discussion, we restricted our analysis 
to the one-dimensional case, i.e., we neglected the cou¬ 
pling between the spin and the motion of the trapped 
bead along the x and y directions on the basis that the 
confinement along such directions is much tighter than 
the trapping along the 2 ) directions. Since the coupling 
terms neglected in Eq. (§ are proportional to 


UJ, 


x,y 


it is worth analyzing the effect of such terms on the dy¬ 
namics predicted in the previous section. In the following 
we will use a perturbative approach, the main scope of 
the treatment being the effect on the interference fringes 
corresponding to the gravity-induced phase difference in 
Eq. (§. In order to do so, we first define the appropriate 
zeroth-order Hamiltonian as 

Hq = D P hjjJxO^a + hjjjyli^h + huJzC^c 

+2AA^(a + a"^) + 2AXy{b + 6’^) + 2AA(c + c^) (12) 

which differs from the Hamiltonian in Eq. 0 for the fact 
that we have added the quantized motion along the direc¬ 
tions X and y, representing by the annihilation operators 
a and 6, respectively. Subsequently the extra gravita¬ 
tional terms due to these transverse oscillation has been 
included through the constants: 


^^x,y — ^x,y ^ 


2moj. 


(13) 


x,y 


where 0^ and Oy are the angles between the direction of 
the gravitational acceleration and the x and y directions 
respectively. 

The evolution of the system with the full Hamiltonian 

H = Ho PV^PVy 

V = {a + at), Vy = \jySy {b + ) , (14) 

can be obtained by treating the terms Vx and Vy per- 
turbatively with respect to the Hamiltonian (12), whose 
eigenstates are products of eigenstates of Sz and dis¬ 
placed number states [26] : 

correspondingly the unperturbed eigenvalues reads: 

E ^ini +Dsl- (16) 


i=x,y,z 


In Eq. (15), we have n = {fix^riy^nz) and D{a) is a 
three-dimensional displacement operator, i.e., D{a) = 
Dx{ax)Dy{ay)Dz{az), with 




-2\sz + 2 AA 2 


1 


2AAa 

HuJt 


2AXy 

huj„ 


(17) 


The first-order correction to the energy due to Vx , Vy 
is zero since 


(s.lWyk.) = 0 . ( 18 ) 


The first nontrivial deviation from the underperturbed 
dynamics can be obtained by correcting the energies up 
to second order and, accordingly, the eigenstates up to 
first order. The corrected energy eigenstates and eigen¬ 
values are therefore given by 



Zp{ 




E 


TT/ 




E, 


(19) 


3^x,y — 


( 11 ) 
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and 


^(1) ^^(0) , y- 

^n,s. ^n,s. ^ (0) _ (0) 

kj^n,s'^ -^n,Sz 


( 20 ) 


where: 


HL 


n,s.-M,s’ = I + Vy) ) , ( 21 ) 


and Z~^ is an appropriate normalization factor. 

For any initial state, the perturbative dynamics at 
any time t can be obtained by expanding it into the 
basis given by Eq. (19) and multiplying each term 


in the expansion by the corresponding exponential fac¬ 
tor exp( —In particular we have calculated 
the evolution of our system from the same initial state 
as in Sec.IL As a result, the motional and spin state 
will be entangled at every time instant t > 0. We 
can quantify how much the perturbed evolution devi¬ 
ates from the unpertutbed one by computing the fidelity 
F = |^^^^(to))| between the perturbed state 

and unperturbed state after one os¬ 

cillation period. We have computed fidelities for a very 
large range of parameters, with A ranging from zero to 0.1 
and 7 a, and 7 ^ ranging from zero to 0.4: this choice in¬ 
cludes also the realistic experimental conditions of our 
proposal, where ujx = ujy = lOuJz^ corresponding to 
7 a, = 7 ^ 0.32. We have found that the fidelity decreases 

with increasing values of either A, 7 a, and 7 ^, but only by 
a very small amount in the range considered: in all cases 
the value of the fidelity stays above 99%, which assures 
that our treatment in Sec. II is accurate enough. This 
is somehow surprising, but the reason relies on the fact 
that the first-order correction in the energies is exactly 
zero due to the selection rules, so that the first significant 
correction to the energies is second order in 7 a, and 7 ^, so 
that even relatively large values of these parameters have 
little effect on the evolution of the spin-diamond system. 

These results remain true when the initial state is ther¬ 
mal. We have indeed found that the fidelity is still very 
close to unity for coherent states whose average thermal 
occupation number is up to 600. This corresponds to ro¬ 
bustness of our schemes up to temperatures of the order 
of ImK, so that only feedback cooling is sufficient for our 
scheme to work El. 


IV. THE EFFECT OF THE RANDOM 
ORIENTATION 

Once the nanodiamond is trapped in the potential, it is 
very likely that the orientation of the NV center quanti¬ 
zation axis, with respect to which the splitting D is com¬ 
puted, is randomly orientated with respect to the trap 
axis. Measuring the optically-detected magnetic reso¬ 
nance spectrum of the NV electron spin in an applied 
magnetic field will reveal the orientation of the NV cen¬ 
ter with respect to the applied magnetic field and hence 


the trap axis. This orientation could be controlled in x-y 
plane by adjusting the linear polarization of the trapping 
light because the nanodiamonds are not spherical. Us¬ 
ing rod-shaped diamonds would increase this control m- 
Alternatively, the birefringence of diamond [28] might be 
used to control the orientation with circularly polarized 
light [29] . 

In this section, without introducing extra measure¬ 
ments and manipulation on the NV’s orientation, we will 
evaluate the systematic errors introduced by the mis¬ 
alignment between the axes of NV center and the trap. 
We will show that the visibility of the interferometry will 
reduce to some degree when these two axes are perpen¬ 
dicular, but in most cases we can get a good resolution of 
the interference fringes. We will essentially use the same 
perturbative methods as in the previous section. 

Assuming we can neglect the coupling with the x and 
y directions, the Hamiltonian we will use is 

H =DSl + huj^c^c + 2AA(c^ + c) 

-2\{c,S, + CySy + c,S,){c^+ c), ( 22 ) 


where , Cy and Cz are the direction cosines of the NV 
symmetry axis z' with respect to the trapping axis z and 
whose values are obtained by a rotation transformation 
between these two axes, bounded by the condition + 
= 1 . Generally there is no exact analytical solution 
for this Hamiltonian, so we divide it into two parts by 
H = Hq i7/, where 

Ho = DSl+nLo,c‘'c+2A\(c‘' +c)-2\c,S^{c^ +c), (23) 
for which we can solve analytically, and 

Hi = -2\{cxSx + CySy){c' + c), (24) 


which will be treated perturbatively. The time evolution 
of any state under Hamiltonian (22) and the effect on the 
Ramsey scheme afterwards can be calculated numerically 
by means of a perturbative expansion to the second or¬ 
der, and based on numerical results we reconstruct the 
fringes of spin zero state population over a range of an¬ 
gles 0 between the trapping axis 2 : and the direction of 
the gravitational acceleration. 

We can investigate the dependence on the orientation 
of the NV center by changing the parameters c^, and Cy. 
Since the system is rotationally invariant along the trap¬ 
ping axis, we simplify the simulation by taking = 0 
and consider c^, varying from 0 (aligned case) to 1 (per¬ 
pendicular case). 

Fig. 3 shows how the interference fringes, given by the 
population Po(^o Ftp) of the state |s ;2 = 0), after the last 
Ramsey pulse has been sent, as a function of the tilting 
angle 0 and of the direction cosine c^, (c^, = 0 corresponds 
to the NV center being parallel to the trapping axis, while 
Cx = I corresponds to the case in which the NV center 
is orthogonal to it). Fig. 3 compares the interference 
fringes in the case of perfect alignment with the fringes 
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in the orthogonal case. The interferometry fringes turn 
out to be qualitatively robust against the misalignment, 
i.e., the visibility of the fringes does not change over quite 
a long range of values of Cx when the NV axis is deviating 
from the paraxial case, and only show a limited reduc¬ 
tion at the extreme case. We can therefore deduce that, 
if we want to perform a proof-of-principle experiment to 
show the existence of the gravity-dependent phase fac¬ 
tors in the evolution of the components of the spin state, 
we do not have to worry too much about correcting the 
orientation of the NV center with respect to the trapping 
axis. 



FIG. 3. Fringes of spin zero population Pq as a function of 
the orientation 0 of the trapping axis, z, with respect to the 
direction of the gravitational acceleration, and of the direction 
cosine Cx = 0 corresponds to the NV center being parallel to 
the trapping axis, while = 1 corresponds to the case in 
which the NV center is orthogonal to it. The initial motional 
state has been taken equal to the vacuum state of the quantum 
oscillator. The other parameters are such that A = 0.01 J and 
AA/cos6> = 11.9 J. 



FIG. 4. Details on how the fringes change when the spin is 
orthogonal to the trapping axis: the visibility reduces but 
does not vanish at the perpendicular case. 


These results in Fig. 3 and Fig. 4 have been obtained 
with A = 0.01 J and AA/ cosO = 10 J, and for a motional 
initial state corresponding to the vacuum state of the 
quantum oscillator. In order to investigate the effect of an 
initial thermal distribution we must carry the simulation 
starting from a generic initial coherent state for the os¬ 
cillator. For computational purposes we have simulated 
the dynamics with the QuTip software m, which as¬ 
sures faster convergence than our perturbative approach, 
and computed the evolution of the system with a generic 
coherent state. Up to an average thermal occupation 
number of CM motion {N) = 600 we have not noted any 
visible changes in the behavior of our system, which as¬ 
sures us that the initial thermal distribution should not 
affect our predictions at least up to temperature of the 
order of ImK. 


V. DISCUSSION AND CONCLUSIVE 
REMARKS 


We will now give the experimental parameters neces¬ 
sary to obtain a good visibility of the interferometric 
fringes in a setup in which we are allowed to vary the 
angle 0. As realistic values for the setup, we consider 
ujz ^ 100 kHz and diamond spheres whose radius R ^ 
100 nm, so that, considering the density 3500 kg/m^ for 
diamond, the corresponding mass is ^ 1.25 x 10“^^ kg. A 
good visibility of interferometry fringes in the population 
in Eq. (10) is given for K = SX AX to/{h‘^ujz cos 0) ^ 10, 
which makes the value of the population change com¬ 
pletely from 0 to 1 when 0 varies between 7r/2 — 7r/20 
and 7r/2 (the ^ axis is horizontal for 0 = 7r/2) (see red 
line in Fig. 3). Assuming that the magnetic field in 
Eq. 0 is generated by a magnetized sphere with radius 
ro = 40 jam and magnetization M= 1.5 x 10^ A/m (typ¬ 
ical for commercial magnets), and zq = 120 /im, we get, 
according to rriz = M • (47r/3)ro and Eqs. ^ and ^ 
we have a magnetic gradient of dB/dz ^ lO^T/m and 
consequently the desired value of K. All these values are 
achievable experimentally and correspond to the values 
of A and AA used to obtain Eigs. 2 and 3. 

With these values, the time necessary to have a com¬ 
plete oscillation of the center of mass is of the order of 
to = 50/is. This must be compared with the typical spin 
dephasing times of NV centers in nanodiamond. NV cen¬ 
ters in isotopically-purified bulk diamond can have elec¬ 
tron spin dephasing times T 2 > 1ms measured with a spin 
echo sequence [31], but such long times have not been 
found for nanodiamonds. The longest times have been 
achieved by making nanodiamonds from pure bulk ma¬ 
terial with a low concentration of spin defects: nitrogen 
impurities must be reduced and ideally also This ma¬ 
terial can then be milled [32| or preferably etched with re¬ 
active ion etching (RIE) to form nanoparticles. Nanodi¬ 
amond pillars with 300-500nm diameter, made with RIE 
have shown a spin-echo T 2 time of over 300 jas [33] . Pillars 
with 50 nm diameter and 150 nm length have achieved 
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a spin-echo T 2 time of 79 fis [34]. By means of appro¬ 
priate decoupling techniques, and by slightly changing 
our Ramsey scheme in order to accumulate the gravita¬ 
tional phase over more than one cycle, such a dephasing 
time can be made larger. With T 2 = 79/is we would get 
perfect visibility of the interference fringes. 

The robustness to the temperatures up to ImK, also 
shows that our proposal is promising for an experimental 
realization. As we have shown, both the terms neglected 
in obtaining our results in Sec. II and the misalignment 
terms are not a problem if we want to perform an exper¬ 
iment aimed at obtaining interference fringes demostrat- 
ing the existence of gravity-induced phases in the spin 
state. 

For other applications of our scheme, such as in metrol¬ 
ogy, we will need a perfect quantitative agreement be¬ 
tween our predictions and the data obtained. In this 
case, we will either have to measure the exact value of 
Cx in order to know to which slice of Fig. 3 our fringes 
correspond, or we will have to correct for misalignment. 
This is currently experimentally achievable, for example 
by the methods shown in Ref. [29] . 

Finally let us stress that our perturbative approach can 
be useful to study the effect of similar spurious coupling 
terms in experimental proposals involving the coupling of 
the motion of nanoparticles with spin degrees of freedom 

miiiiiiss]. 

In this paper we have studied the dynamics of an 
NV center in a harmonically trapped nanodiamond, 
considering the effect of unwanted coupling between the 
spin and the directions orthogonal to the direction where 
the conditional displacement takes place, inducing a 
gravity-dependent phase which can be detected by spin 
measurements only. We have shown that a perturbative 


treatment of the additional terms in the Hamiltonian 
proves that, for the experimental parameters charac¬ 
terizing the setup proposed, the fidelity between the 
realistic and the unperturbed state is always above 99%. 
This assures that the experiment should confirm the 
one-dimensional theory, which allows for a very clear 
physical interpretation of the results. Moreover, also the 
robustness with respect to thermal fluctuations in the 
initial state is retained, which is the most important 
feature of our proposal, since no ground state cooling 
is required, only some feedback cooling to mK temper¬ 
atures in order to guarantee that no anharmonic effects 
come into play in the dynamics. The same perturbative 
approach can be applied to treat the misalignment 
between the NV center and the trapping axis. In the 
case of misalignment, the quantitative agreement with 
the one-dimensional model is lost, but the qualitative 
behavior of the interference fringes stays valid up to 
very high deviations from the perfectly aligned case. 
Therefore, a proof-of-principle experiment aimed at 
showing the existence of Ramsey interference fringes 
induced by gravity can be performed without caring too 
much about the alignment between the NV center and 
the trapping axis. 
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